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ABSTRACT 


Both analytical and experimental studies were conducted on a pro- 
pagating crack in a viscoelastic material. By applying an extension of 
the "correspondence principle’ the stress and displacement at the crack 
tip were found as functions of the crack tip stress intensity factor. 
An energy balance was made for the specimen which allowed this factor 
to be related to geometric and load parameters. Fracture characteriza- 
tion was then performed by experimentally relating the crack tip stress 
intensity factor to the crack velocity. 

The theory was applied to solid propellant fracture tests; however, 
the dewetting behavior of the propellant was not accounted for by the 


theory. 
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I. INTRODUCTION 


In the elastic crack propagation problem, it is useful to define a 
crack tip stress intensity factor, K. This factor, which includes the 
effects of specimen loading and geometry, gives a measure of the stress 
concentration at the crack tip. As the load on a Specimen containing a 
stationary crack is increased, K increases until a critical value is 
reached, whereupon the crack rapidly propagates. This critical value, 
Ke» 1S experimentally measured and found to be a material property. | 
A material may then be characterized by this constant, for if any com- 
bination of geometry and loading should produce Kae crack propagation 
will occur. 

For cracks in viscoelastic materials, Knauss and Mueller [Ref. 2] 
showed that in effect there 1s a certain minimum loading which will 
initiate crack propagation; that is, there is a loading below which the 
crack velocity is small enough to be considered equal to zero. However, 
the velocity of propagation will vary over a range of several orders of 
Magnitude, increasing with increasing load. Thus the auestion to be 
answered in the viscoelastic fracture problem is not simply will the 
crack propagate, but rather at what velocity will it propagate. Knauss 
[Ref. 3] proposed that for a given viscoelastic material, K(¢) is a 
unique function, where ¢c refers to the crack velocity. Therefore it is 
natural to experimentally characterize viscoelastic materials by this 
relation. 

"The material property, K-,.iS unique for most materials and geo- 
metries. However, there iS a variation of K_ with initial crack length 


(usually small enough to be ignored) explained by consideration of the 
growth resistance curve as shown in Ref. 1. 





A common method of obtaining viscoelastic solutions is to use the 
"correspondence principle.” Lee [Ref. 4] showed that viscoelastic 
boundary value problems possess the property that an appropriate elastic 
solution in the Laplace transform domain may be inverted into the time 
domain to give the viscoelastic solution. However, .this "correspondence 
principle’ does not in general apply to problems where the boundary 
conditians are mixed and apply over regions which change with time. 

Now looking at the problem of the propagating viscoelastic crack, 
note that the crack edge boundary condition immediately behind the crack 
tip 1s zero stress while the boundary condition on the extended line of 
propagation of the crack, immediately in front of the tip, is zero 
displacement. Propagation of the crack implies that portions of the 
boundary under zero displacement are continually changing into zero 
stress boundaries. Thus, the restriction on the ‘correspondence 
principle" mentioned above prohibits its application to this problem. 

A different solution method has been used by Dietmann, Knauss and 
Mueller [Refs.2,5,6,7] in their various studies of propagating cracks 
in viscoelastic materials. In all of these cases, the actual problem 
was approximated by one in which crack propagation was a step-wise 
process. The resulting solutions in all cases contained the length 
of the crack increment used in the step-wise process as a material 
property; however, this length does not physically exist in the actual 
problem and therefore seems to be only the result of the solution 
technique. 

In the solution of the viscoelastic fracture problem developed in 
this paper, the "extended correspondence principle" will be used. This 


extension is due to Graham [Ref. 8] who extended the classic 








“correspondence principle" to cover several cases where mixed boundary 
conditions applied over regions which changed with time. Using this 
extension, both the viscoelastic stress and displacement in the 
vicinity of the tip of a continuously propagating crack will be 
obtained as functions of the stress intensity factor, K. The remainder 
of the salution will consist of applying an energy balance tothe 
specimen and relating K to the loading and geometric parameters of the 


specimen. 








ie tmeORY 


Consider a thin, narrow sheet of isotropic, homogeneous, linear 
viscoelastic material loaded uniaxially by o,(X,t) (where X refers to 
geometric coordinates and t to time) containing a continuously pro- 


pagating center crack of length 2c(t) as shown in Figure 1. 


y 
q(X,t) 


is 


Figure 1. Specimen Geometry and Loading. 


Assuming that the material is thermorheologically simple, the analy- 
sis for the stress intensity factor is made at a constant reference 
temperature, To: This solution will then be valid at any other 


temperature when appropriately shifted. 


A. CRACK TIP STRESS AND DISPLACEMENT 

In this analysis a knowledge of the work produced at the crack tip 
will be necessary. This further requires an evaluation of the stress 
and displacement fields in the vicinity of the moving crack tip, which 


1S now undertaken. 








If the sheet and loading are doubly symmetric, a single quadrant 
may be considered, as shown in Figure 2, where the lower right quadrant 
is removed and replaced by the stresses which act along the sheet's 


horizontal centerline. 
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Figure 2. Single Quadrant of the Doubly Symmetric Specimen. 
Let the entire boundary of the sheet be defined as 8 with Bo (t) 
representing the uncracked portion of the x-axis and 8, (t) representing 


the remaining boundary. Then 8 = 8, (t) UJ 85(t). The boundary conditions 


are: 
o.(X,t) = Q on B 
0 on crack edge 
o(X,t) =< Q on x = 0,1 on 6, (t) (1) 


o,(X,t) on y=b 
u(x,t) = 0 on Bo (t) 
where u refers to displacement and the subscripts s and n respectively 


refer to normal and shear components. The boundary conditions of 








Equation (1) meet the criteria of Graham's "extended correspondence 
principle." 

However, to apply this correspondence, the elastic stresses and 
displacements in the vicinity of a stationary crack tip are needed. 


They have been found by Irwin [Ref. 9] to be: 





1 G ch 36 
cos 5 | + sin a Sin , 2 (2) 
2 





aap ia 8 7 - Ce) 
Un= Go oe =Sin 5 2 2v - cos 7 (3) 


where the subscript 1] on the stress intensity factor indicates that the 
crack 1S propagating in the opening mode as defined by Ref. 10. 

Setting 6 = 0 in Equation (2) and 6 = 7 in Equation (3) gives the 
stress and displacement, close to the tip, on 8, (t) and 8, (t) 
respectively. It is reasonable to assume that Equations (2) and (3) 
apply to slowly propagating cracks where the inertia effects are neglig- 
ibly small in the equilibrium equation. Changing to an inertial 
coordinate system, where x = c(t) + x' (refer to Figure 2) and assuming 
the incompessible case (vy = 5) the elastic stress and displacement in 


the vicinity of a slowly moving crack tip become the following: 


K,(t) ; | 
o® (X,t) = ———_ h (¢ - t) on B(t) (4) 
¥ 2n (x-c(t)) 
eel) 
uP (Xe) = 425 c(t)-x h(t-t) = 
edge 
us OG) on x = 0,1 8, (t) (5) 
Us (X,t) | on y = b 


Cth elastic theory, the stress intensity factor in these equations 
is the critical one, Kic, since the crack is moving. However, in the 
viscoelastic case, with continuous propagation, this definition of the 
critical value loses its significance. 
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where superscript e refers to elastic solutions, us (X,t) and UD (Xa) are 
elastic displacements which are functions of geometry and modulus, h(t) 
is the unit step function and t is the time at which the crack tip 
passed a point on the horizontal axis located at x. 

Graham has shown that Equations (4) and (5) are the solutions to the 
viscoelastic problem in the Laplace transform plane if the elastic 


modulus is replaced by the transform modulus, E>pE..(p). Therefore, 


rel 
the field quantities of interest are: 


- Kealet) . = 
* | 
oF (X,p) = p= h (t - 4) on Bo (t) (6) 
* ae ieee | 
Uy (X;p) = 2 z= pea] D ie c(t) - X h(t = | Cy) 


- ; on 8, (t) 
where: foe Pir(tydt = [f(t)] 


0 


Since pD..(P) = 1/pE...,(p), where D., refers to creep compliance, 


rel 
the viscoelastic solution inverts back into the time domain with the 
stress distribution identical to the elastic case and the displacement 


at the crack tip edge given by: 





U, (x,t) = 24 é Bt K, (0) c(0) - x h(0-t) 


Since h(0 --t) = 0 and h(t - t) = 0 fort < t, Equation (8) may be 


rewritten as: 


t 
Uy (xat) 2 7 De-(t -r) aie c(t) - «|e (9) 


1] 





Time t is seen to be equal to time t minus the time which it took 
for the crack to run between x and c(t), that is: 


a C12) dee (10) 


Cm 


ct > 
i 

iT] 

tt 
ctr 


where C, is the mean velocity in the interval [x,c(t)]. Thus, time t 
accounts for the history which has.affected the displacement at point x. 
Or to put it another way, the crack displacement at point x, at the 
present time t, has been creeping open ever since the crack tip passed 
that point at time t. 

The expressions for the stress and displacement, which have been 
found in this section,will allow the energy lost at the moving crack 
tip to be evaluated. Thus, the analysis now proceeds to an energy 


balance. 


B. ENERGY BALANCE 

An energy balance will now be applied to the entire specimen. First 
the specimen energies are defined as follows: W is the work input to 
the specimen by forces acting on the boundaries; I is the reversibly 
stored energy in the viscoelastic specimen; D is the energy which 1s 
dissipated throughout the specimen due to the viscous action of the 
material; finally, S is the energy which is dissipated at the crack 
tip in the creation of new surface area. 

Neglecting heat and inertia terms, the First Law of Thermodynamics 
may be written: 


Ws T+ D +S 


Z 





The significant variables in the problem are c and t, thus the power 


equation becomes: 


























aw aw | _ al a} , 2a oD | , asi. 2S 
meio oe | ~~ sae) © * oe) * ac | © ae] Tc | © tt 
c C 
Rearranging terms the following equation is obtained. 
am 6ot «|. aD | as - | Cities Oded Dis inc-0 3) 
-# 3¢ | _ Pea S| 1) 
ie C 


The left hand side of Equation (11) may be thought of as describing 
the elastic behavior of the viscoelastic speciman (that behavior which 
is independent of time) while the right hand side describes the 

aS z : 
a arie O since there 
1S nO mechanism to create energy S. Thus, for a fixed crack length, 


viscoelastic (time dependent) behavior. However 


the first law becomes W = I + D and the right hand brackets of 
Equation (11) sums to zero. Also, viscous dissipation is physically 
an energy which must occur over finite time, thus in the left hand side 


of Equation (11): 





Therefore if time is held constant, Equation(11]) becomes: 
ol oI J9 : 
Fg ac ae | : Me) 
1 


Of the energies which appear in Equation (12), W is the most 
directly evaluated. It is simply the energy produced by the applied 
force, O. moving through the displacement at the boundary. The 


remaining energies of Equation (12) will now be considered in detail. 


1G) 





ee neversiuly stored Energy 


To further clarify the definition of I at any time t, consider 
the following. Uniaxial viscoelastic behavior can be modeled by the 
three-element mechanical model consisting of a Maxwell element and a 
spring loaded in parallel as shown in Figure 3. The short-time spring 
constant for the model is EY Sob t E. where E >> E. The spring 
constants of the model correspond to the slopes in a stress-strain 
plot while model force and displacement correspond respectively to 


viscoelastic stress and strain. 





Figure 3. Stress-Strain Plot and Mechanical Model for Visco- 
elastic, Reversibly Stored Energy. 


Consider the loaded viscoelastic body (model) at any time t 
(point A). Allow an instantaneous relaxation, thus losing no energy 


in viscous dissipation in the model dashpot, and calculate the energy 
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released. Since this is a conservative process, the energy released 


will be the reversibly stored energy. 
The relaxation occurs with the short-time modulus to a point 


af zero stress (point B). The energy releasedper unit voiume is the 
area under the path AB in the stress-strain plot or: 
; 2 
mie es onlt) 


VOL Cee 
g 


Comparing this to the strain energy of a uniaxially loaded elastic 
body: 


aU yen 
VOL 2 


the following correspondence may be made between elastic strain energy 
and internal energy in specimens of identical geometry which are loaded 


uniaxially. 


U+ ] 
if Eset (13) 
' g 
and o + (t) 


cm EWerGy Dissipated at thes Crack Tip 
The* change in energy S with respect to length, at the crack tip, 


required in equation (12), is now evaluated. Holding t = t., allow a 


Go 
small propagation of a (see Figure 4). Since this is for constant 
time, the behavior of the viscoelastic material 1s elastic and the 
evaluation technique used by Irwin [Ref. 9] for elastic materials may 


be applied. Thus, the energy lost in this process is: 











Figure 4. Crack Tip for Incremental Propagation, a, at t= Tt. 


Substituting for the maximum values of displacement and stress, this 


becomes: 


—s 


a 
AS = ae UL ‘Coma te) On Gee te) dx' 


or in inertial coordinates: 


Equation (14) requires that the displacement of Equation (9) be 
translated by x = a. Note that the lower limit on the integral of 


Equation(10), which is a function of position, thus changes to: 


t +22 fires t- A 
m m 


(15) 
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With this change in limits, the displacement, Equation (9), and the 


Stress, Equation (4), are substituted into Equation (14) to give: 
t_)+o i 
K, (t_) c( € 2 
oa f f Derltert) og 
= OT 
C 


ts vc(t) - x ms drdx (16) 





Consider the compliance in Equation (16). At the upper limit 

in time it becomes D.-(9) while at the lower limit it is: 
[ - X + - 
Der 
Ch 

which 1s a function of x. Now evaluating at the upper limit in x, it 
again is D..(0) while at the lower limit is becomes Dap lo/e,,) « Thus 
since a is small, the compliance may be approximated by the constant 


value, D..(U). The integrand in time is now a perfect derivative and 


or! 


Equation (16) becomes: 





oc = Y = 
y= TE po) |i) TET =e 


- K, (t - fal =) yclt - -x + ’ dx (17) 


~ 


Consider the integrand of Equation (17) which is evaluated at 
. ; A“ Q : 
any x for which c(t) <X< c(t.) +a. Choosing x, c(t - =) is seen 
m 
to be equal to the crack length at the time it passed x minus the length 


of time it took to propagate through a distance o at its. mean velocity, 


Co In other words, c(t - a is equal to x - a and the second term in 
m 





the brackets vanishes. Now the time dependence in this energy has been 
evaluated and Equation (17) may be conveniently written in local 
coordinates as: 


2 O 
i, (t_) : 
oo = § a [ gp 
cr x! 
0 


je 





which when integrated becomes: 





Since a is small, the change in S with respect to c is: 


2 
as _ as _ Kyit) DCO) 


J ne (18) 


where the subscript on t has been dropped since t was a typical time. 
2. Specitic Specimen Geometry 
At this point, the energies required in Equation (12) have been 
defined or evaluated. However, since W and I depend on the particular 
geometry and loading, the evaluation for these energies has been only 
general. A specific case is shown in Figure 5 where the applied load, 


o(X,t), is: 


a(x, thet) = gy tO ee iy, (19) 


Oe formmwet) 

The strains in the specimen may be thought of as a superposition 
of the uniaxial strain for x 2 c(t) and the strain produced by the 
stress field in the vicinity of the crack tip. Since the uniaxial 
displacement of the boundary edge is e(t)b, choosing b sufficiently 
large will make this displacement of higher order than that produced on 


the boundary by the stress field at the crack tip. Thus holding 


18 
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Figure 5. Specimen Geometry and Loading. 


time constant and allowing an incremental propagation of dc, while 
removing load from the boundary as specified by Equation (19), produces 
negligible displacement, hence negligible work at the boundary. 
Therefore: 


oW 


Na = 0 (20) 





5 

For this specimen, if the uncracked length 1 - c(t) is maintained 
large enough so that a strip at the right hand edge is under uniaxial 
stress, that is, unaffected by the crack tip stress field, Ref. 1] gives 


the change in the elastic strain energy as: 





2 
aU 4. eOue 
vs =~ b (0) 
t 
Thus by applying Equation (13): 
; 2 
1 | 
ae = - Fb (22) 
t g 
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Substituting Equations (18), (20) and (22) into Equation (12), noting 
that now surface energy is released from two surfaces and that ES = 
m7 De AO); gives: 

K, = o. Yb (23) 


Note that since o, is constant, Ky 186 a constant. This result 


1s identical to that which is obtained in the elastic case. 


C. ADDITIONAL RESULTS 

Lindsey [Ref. 11], applying Graham's “extended correspondence 
principle, transformed the elastic expressions for stress intensity 
factors directly into the viscoelastic case. His results agree with 
Equation (23) and include the additional case of a strip with 
displacement boundary conditions under a constant displacement, u, 
(see Figure 6). (Uncracked specimens of this geometry are commonly 


used in the propellant industry and are referred to as “biaxial strips" 


Since they are in a state of biaxial stress.) Lindsey's solution is: 


(t) yr (24) 


a(t) 


SR RRR eeees: 


RM ©AMhO~AwAs ————— 





S| _—e a 


Figure 6. Cracked Biaxial Strip. 
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D. TIME-TEMPERATURE SHIFT 

The expressions for K,(t), Equations (23) and (24) are for a 
constant reference temperature, ae Since the viscoelastic material 
was assumed to be thermorheologically simple, Ref. 12 shows that the 
modulus, hence the stress, may be shifted from any temperature, T, 
to the value at ih by multiplying it by the factor lip Elie where 
o 1S the material density. The time scale 1s corrected by dividing it 
by a time-temperature shift factor, ay SInee sp] Po? the expression 


for the stress intensity factor become: 


7 ir t 
K(E) - Oo. f= yp SUVeSS. Bete (25) 


Baie Oe (Pah. 
K, (=) = = “( ==) q3 b Bea. tee. (26) 


These expressions can now be used to guide experimental measurements 
of kK, (t/ar). For, in light of the previous discussion concerning the 
relationship between Ky and ¢, laboratory measurements of these two 
parameters in fracture tests will constitute a viscoelastic fracture 


characterization of the material. 
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ITI. EXPERIMENTAL RESULTS 


Experimental testing was conducted on a CTPB solid propellant. 
This filled material consisted of solids (oxidizer) comprising 
approximately 85% of the total volume, bonded in a rubbery (fuel) 


matrix. 


A. CONSTANT STRESS CASE 

The specimen configuration to achieve the stress boundary conditions 
of Figure 5 is shown in Figure 7. The specimen was cut from samples 
designed for use in biaxial tests with a gage length of 1.0 inch and 
thickness of 0.1 inch. The wooden tabs were segmented with saw cuts 


which stopped in the adhesive bond just short of the propellant. 


lA, ~ 
i | Redwood 
© © wt : ee tabs 
W 
d | 
Adhesive — : : 
bond 
Propellant 


c(t) ———>» 


© 2 G ° O p D 
’ 


(Full scale) 


WWWWw WW W W 
Figure 7. Specimen Configuration for Stress Boundary Conditions. 


Ze 





(To penetrate through the bond into the propellant would have created 
a stress riser and a point of probable crack initiation.) The width, 
w, Ideally would be very small, allowing shear-free edges; however, 
practical considerations forced a-width of w =-0.25 inch. 

Initial experiments, with the right-hand edge of the propellant 
flush with the wooden tab edge, showed that soon after loading a fast 
running crack would initiate at a local defect along the right edge of 
even a carefully cut specimen. Adding a 1/8 inch unloaded length reduced 
the stress along this. edge and eliminated the problem. To compensate 
for this departure from uniform stress at the right hand edge and yet 
stil] meet the restriction leading to Equation (21), Stee was limited 
to 2.0 inches. 

The thickness of the saw cuts inthe wooden tabs was large enough to 


al low unrestrained Ey» or d/w > e The upper segments were 


X max’ 
individually suspended and the lower segments were individually loaded 
by weights W. An initial crack of one inch was cut with a razor blade. 
A specimen was tested by applying all loads at t = 0. The length 
of the propagating crack, c(t), was measured with an optical comparator. 
As the crack passed a segment, the weight on that segment was removed 
in proportion to the width of the segment traversed by the crack. The 
tests were conducted in an environmental chamber with the relative 
humidity less than 30%, at six temperatures ranging from 15°F to 85°F. 
The stress intensity factor for each test was calculated by using 
Equation (25), assuming op = 1, and plotted on a graph of stress 
intensity factor versus crack velocity, for each specific testing 
temperature (see Figures 8 through 13). In this way the shift due to 


temperature, in the ordinate, was accounted for. However, complete 


2S 





shifting of the function required that ay be known. This was obtained 
graphically by shifting the curves for various temperatures along the 
abscissa, relative to the curve at Ty» so that a continuous function 
was obtained. The resulting characterization curve is shown in Figure 


14, where i = 7@LF . 


B.. CONSTANT DISPLACEMENT CASE 

Additional tests were conducted on cracked biaxial strip specimens. 
Though the ends were free, rather than restrained as shown in Figure 6, 
the long length relative to the width allowed the actual specimen to 
closely approximate the theoretical one. The strips were six inches 
in length, one inch in width (2b) and 1/10 inch thick. Redwood tabs 
bonded to the horizontal edges allowed for gripping. An initial crack 
of 0.75 inch was cut, as suggested by Mueller [Ref. 7]. As béfore, 
the length of the propagating crack was measured as a function of time 
with an optical comparator. 

The testing was performed with the apparatus which was used for 
Similar tests on solithane, as described in Ref. 13. Briefly this 
consisted of an Instron Universal Testing machine equiped with special 
jaws to clamp the specimen and apply displacement u. while maintaining 
the upper and lower boundaries parallel. The displacement was applied 
at rates ranging from one inch/minute for low temperatures to five 
inches/minute at the high temperatures. A conditioning chamber al lowed 
the tests to be conducted at the same specific temperatures as in the 
previous case. 

Using Equation (26), K, was calculated and plotted in Figures 8 
through 13. Shifting was performed as in the previous case, with ie ~ 


70°F, and the master characterization curve plotted in Figure 14. 
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IV.. DISCUSSION OF RESULTS 


In the stress boundary condition case, the crack tip stress 
intensity is constant. This constant stress intensity should produce 
@ constant velocity crack, a fact which was verified in these experi- 
ments. Therefore, each sample tested, for a given loading and 
temperature, produced a unique point on a plot of “K,(¢). 

However, for the displacement boundary condition case, stnce the 
Stress relaxes, the crack will decelerate until a (rubbery modulus ) 
is reached. For the propellant tested, this relaxation time was on the 
order of a year. Therefore, for all experiments of this case, the 
crack was constantly decelerating. Thus Ky (¢) was a continuous function 
and in this study was shifted more readily than the stress data, since 
for the later case, limited samples precluded enough testing to obtain 
a well defined function for each temperature. For this reason the 
Shift factors obtained in the displacement case were used in the stress 
case to shift the data to the master curve shown in Figure 14. 

The theory of Section II considered a purely viscoelastic material. 
However, the solid propellant was a filled material, which when loaded, 
exhibits the additional mechanism of dewetting. (Dewetting is the 
breaking of bonds between solid filler particles and the rubbery 
matrix.) Wood [Ref. 14] has shown that propellant dewetting is not a 
rate process, rather it occurs upon initial application of stress. 

This behavior has a significant bearing on the experimental results 
obtained in this study as will now be shown. 

In both cases tested, the measured load parameter was stress, thus 


automatically accounting for the modulus relaxation in the displacement 
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case. However, the two stress histories were quite different as shown 
in Figure 15. In the stress case, stress was constant while in the 
displacement case there was an initial extremely high stress which 


continual ly relaxed.> 


tg 


Constant Stress Constant Displacement 


Figure 15. Stress Histories. 


Thus a sample tested in the stress case exhibited minimum dewetting 
because of its history of constant stress. However, in the displacement 
case, for a corresponding stress at t 2 10t. (where t. is the time 
required to apply the load) much more dewetting had occurred due to the 
initial high stress. Consequently the crack propagated through a 
material which had been significantly altered from that in the stress 
case. The propagating crack had fewer bonds to break per unit area, 


compared with the stress case, and accordingly its velocity was as much 


3The theory of this paper does not account for any transients, 
rather considers only step loading. However, the application of the 
loading in the displacement case required finite time. Therefore, data 
taken for this case were taken at t 2 10tqg assuming that at this time, 
any effects of the loading method had died out. 
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as two orders of magnitude higher. This produced significantly 
different characterization curves for the two cases as shown in 
Figure 14. 

Since dewetting produces microscopic voids in the material, stiff- 
ness is reduced as dewetting increases. Using this fact, the increased 
dewetting of the displacement case was verified in preliminary tests at 
15°F. In these tests, for identical stress and time, the modulus in 
the displacement case was one-half that of the stress case. 

In an attempt to reduce the initial high stress peak in the dis- 
placement case, thus reducing the dewetting, additional preliminary 
tests at 15°F were conducted. The maximum stress was limited to that 
which produced stress intensity factors comparable with those of 
the stress case for the same temperature. Equivalent moduli showed 


that the amount of dewetting was in fact equivalent. However, as soon 


Gone 


Thus even if the crack propagated at the same velocity as in the stress 


as the displacement rate was stopped, the stress relaxed due to E el 
case (approximately 0.1 inch per hour) relaxation of the driving 
Stress decelerated the crack to zero velocity over distances too 
small to be accurately measured. Thus, no additional characterization 
data were obtainable with this loading method. 

The accuracy of the measured velocity in the stress case was con- 
Sidered superior to that of the displacement case, especially for low 
velocities. For in the displacement case, c(t) had to be measured over 
relatively short periods of time to obtain an accurate crack propagation 
function. (Velocity was obtained by graphically measuring the slope of 


this function.) As ¢(t) became small, Ac(t) for the time interval 
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mentioned above also became small. Thus, because the experimental crack 
propagation was irregular due to inhomogeneities in the propellant, the 
actual velocity was obscured as those irregularities approached the 
magnitude of Ac{t). However, in the stress case, ¢(t) was a constant 
and. Ac{t) could be measured over large time intervals. 

The scatter in the characterization curves is due to the inaccur- 
acties mentioned above for the displacement case and also to the 
imhomogeneous nature of a highly filled viscoelastic material such as 
propellant. By comparison, similar characterization of solithane, an 
unfilled and relatively more homogeneous material, by Francis et. al. 
[Ref. 13] showed significantly less scatter.” 

Referring to the characterization curves, in both cases there appears 
to be a K critical below which a crack will not propagate. Again for 
both cases, especially apparent in the stress case, there appears to be 
a: maximum velocity of propagation; however these maximum velocities 
are somewhat misleading. Preliminary tests showed that at very high 
loads the failure was general and catastropic, initiated at voids and 


other stress risers throughout the specimen as well as at the crack tip. 


+ Solithane at room temperature has a relaxation time of less than 
one minute. Thus, though tested for the displacement case, essential ly 
constant velocity crack propagation in these tests allowed the additional 
advantage of increased accuracy of the measured crack propagation 
function for reasons mentioned above. 
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V. CONCLUSIONS 


The theory presented in this paper permits’ fracture characterization 
of a viscoelastic material to be made. However,.for a filled visco- - 
elastic material such as a solid propellant, dewetting behavior affects 
the crack velocity. (The stress boundary condition case is the limit 
case of minimum dewetting.) Additional theory is necessary to ACeourt 
for dewetting. This theory would then be expected to provide for a 
superposition of the displacement characterization curve upon the stress 


characterization curve (See Figure 14). 
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